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Abstract 

Using experiments and simulations, we investigate the clusters that form when colloidal spheres 
stick irreversibly to - or "park" on - smaller spheres. We use either oppositely charged particles 
or particles labeled with complementary DNA sequences, and we vary the ratio a of large to 
small sphere radii. Once bound, the large spheres cannot rearrange, and thus the clusters do 
not form dense or symmetric packings. Nevertheless, this stochastic aggregation process yields a 
remarkably narrow distribution of clusters with nearly 90% tetrahedra at a = 2.45. The high yield 
of tetrahedra, which reaches 100% in simulations at a = 2.41, arises not simply because of packing 
constraints, but also because of the existence of a long-time lower bound that we call the "minimum 
parking" number. We derive this lower bound from solutions to the classic mathematical problem 
of spherical covering, and we show that there is a critical size ratio a c = (1 + y/2) « 2.41, close 
to the observed point of maximum yield, where the lower bound equals the upper bound set by 
packing constraints. The emergence of a critical value in a random aggregation process offers a 
robust method to assemble uniform clusters for a variety of applications, including metamaterials. 



Understanding the geometry of clusters formed from small particles is a fundamental 
problem in condensed matter physics, with implications for phenomena ranging from nu- 
cleation pQ to self-assembly [2]. Colloidal particles are a useful experimental system for 
studying cluster geometry and its relation to phase behavior [3j for several reasons: they 
are large enough to be directly observed using optical microscopy; their assembly can be 
understood in terms of geometry [HE]; and they can be driven to cluster by a variety of con- 
trollable interactions, including capillary forces [2J, depletion [6], fluctuation- induced forces 
[7] , or DNA- mediated attraction [B] . Colloidal clusters are also useful materials in their own 
right. They can be used, for example, as building blocks for isotropic optical metamaterials 
known as metafluids [9HTT]. Tetrahedral clusters are of particular interest for metafluids 
since the tetrahedron is the simplest cluster with isotropic dipolar symmetry [9]. An un- 
solved challenge for this application is to determine the interactions and conditions that 
enable assembly of bulk quantities of highly symmetric, uniform clusters such as tetrahedra. 

With this motivation in mind, we study experimentally the geometry and size distribu- 
tion of binary clusters formed when small colloidal spheres are mixed with an excess of large 
spheres that stick irreversibly and randomly to their surfaces (Figure [T^i). An obvious way 
to control the cluster geometry in such binary systems is to vary the size ratio. One might 
expect that at certain ratios the particles could arrange into dense clusters or "spherical 
packings" - arrangements of spheres around a central sphere that maximize surface density 
P2HII]- Such packings have long been used in modeling the microstructure of dense, dis- 
ordered atomic systems [151 HE] ■ But unlike atoms, colloidal particles can stick irreversibly, 
such that two particles bound to a third show no motion relative to one another. This type of 
binding occurs frequently in strongly interacting, monodisperse colloidal suspensions, which 
consequently form fractal aggregates instead of dense glasses [T7J HH] . Similarly, in the binary 
systems we study, the irreversible and stochastic process of sticking precludes the formation 
of dense or symmetric packings. The large spheres park, rather than pack, on the surfaces 
of the small spheres. 

Surprisingly, this random and non-equilibrium process can produce clusters of uniform 
size. Our experiments show that at a size ratio a = Rbig/Rsmai\ = 2.45, where i?bi g and 
-Rsmaii are the sphere radii, nearly all of the clusters contain four large spheres stuck to a 
smaller sphere (Table [T]). In these experiments we use a 100:1 stoichiometric ratio of the 
two sphere species, statistically ensuring that each cluster contains only one small sphere 
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FIG. 1. (a) Two colloidal sphere species are mixed together to form clusters, (b) Oppositely 
charged polystyrene spheres cluster due to electrostatic attraction. Optical micrograph shows a 
tetramer (iV = 4). (c) Polystyrene spheres labeled with complementary DNA strands (not to scale) 
cluster due to DNA hybridization. Optical micrograph shows a trimer (N = 3); the small, central 
sphere is fluorescent. 



surrounded by two or more larger spheres. After waiting several days for the average cluster 
size to saturate, we measure the distribution of N, the number of large spheres bound to 
each small sphere [19J. We do not count single large spheres, nonspecifically aggregated 
clusters of large spheres, or clusters with multiple small spheres. While there are many 
isolated large spheres due to the high stoichiometric ratio, the latter two types of cluster are 
rare. 

The A^ = 4 tetramers that we observe are not dense packings or, in general, symmetric 
arrangements. As can be seen from the images in Figure [TJ there is space between the large 
particles, and the resulting tetrahedra are irregular. Moreover, the ratio a = 2.45 is well 
below the value a = 4.44 found by Miracle et al. [20] for efficient tetrahedral packing in 
binary atomic clusters. In fact, at a = 4.29, closer to this bound, we see much smaller 
clusters and few tetrahedra. The sparsity of large spheres in the clusters is a result of the 
irreversible, non-equilibrium, random binding: once the big particles stick to the smaller 
ones, we do not see them detach or move relative to one another. We expected such a 
stochastic process to lead to a much broader distribution of clusters. At other values of a it 
does (Table [T]), but at a = 2 Ah we obtain 90% tetramers. 
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Size ratio a 


1.94 2.45 3.06 4.29 


N = 6 


6.3 0.0 0.0 0.0 


N = 5 


39.2 0.8 0.0 0.0 


N = 4 


54.4 90.2 18.6 0.7 


iV = 3 


0.0 6.6 69.9 35.9 


N = 2 


0.0 0.8 10.9 51.0 


iV = 1 


0.0 0.8 0.6 11.1 


iV = 


0.0 0.8 0.0 1.3 



TABLE I. Experimentally observed cluster size distributions for charged colloids. Percentages of 
total are listed. The distribution for a = 2.45 (red) is sharply peaked at N = 4. 

The high yield of tetramers occurs in two experimental systems with different types of 
interactions. In both systems the interactions are specific, strong, and short-ranged, and 
the particles do not rearrange once bound. In the first system the clustering is driven by 
electrostatic interactions. We mix large, positively-charged particles with small, negatively- 
charged particles, as shown in Figure [T]d. To adjust a, we use several different particle 
sizes [19]. We add salt to reduce the Debye length to approximately 3 nm, small enough 
to ensure that the interaction range does not significantly influence the effective particle 
size. In the second system the clustering is driven by hybridization of grafted DNA strands 
|19j . As shown in Figure we mix small and large spheres labeled with complementary 
DNA oligonucleotides [21]. We work well below the DNA melting temperature so that the 
attractive interaction is many times the thermal energy [22J. 

To better understand why the distribution is sharply-peaked at = 4 for a = 2.45, we use 
simulations and analytical techniques that account for the irreversibility of the aggregation 
process. Our simulations use a "random parking" algorithm [23H26] to model the formation 
of clusters. The algorithm involves attaching large spheres to randomly selected positions on 
the surface of a small sphere, subject to a no-overlap constraint [19J. We do not model the 
finite range of the interactions, which in both experimental systems is small compared to the 
particle size, or the diffusion of the particles prior to binding. In accord with experimental 
observations, the particles are not allowed to rearrange once bound. We repeat the process 
numerically to obtain distributions of cluster sizes as a function of a single parameter, a. 
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FIG. 2. (a) Yield curves, as determined by simulations, for iV-particle clusters, 2 < N < 8, where 
the critical size ratio a c is marked with a black line. Below are histograms for (b) DNA-labeled 
particles (blue) at a = 1.90 and (c) charged particles (red) at a = 2.45, as observed in experiments 
and as predicted from simulations (gray). Error bars are 95% confidence intervals (Wilson score 
interval method). 



The simulations find a 100% yield of tetramers at the size ratio a ~ 2.41. As in the 
experiments, the large particles in these tetramers are not densely packed, and the clusters 
are therefore distorted tetrahedra. We also find that while the yield of any particular cluster 
can be maximized by varying a (Figure [2^l), the yield approaches 100% only for dimers 
(N = 2) and tetramers (N = 4). Interestingly, the yield curve for tetramers has a cusp at 
its peak, showing that the size ratio a c at the maximum is a mathematical critical point. 

The simulated distributions agree well with those found experimentally (Figure [2]d,c) 
for both electrostatic and DNA-mediated interactions. For instance, at a = 2.45 with 
electrostatic interactions, we find a sharply-peaked distribution consisting almost entirely 
of tetramers. This value of a is close to but not precisely at the critical value, so a small 



yield of trimers is predicted and observed experimentally. In contrast, at a = 1.90 we find 
a mixture of mostly N = 4 and N = 5 clusters in both the DNA system and simulations. 
Some discrepancy arises between the simulated and experimental histograms because the 
yield curves in Figure [2^, are steep; a slight error in the effective size ratio can shift the 
cluster distribution. Nevertheless, the random sphere parking model successfully reproduces 
both the large yield of tetrahedra near a c and the details of the measured histograms at 
various other a. 

That we can reproduce the same phenomenon in two different experimental systems and 
in a one-parameter model suggests that the critical size ratio a c has a universal, geometrical 
origin. Intuitively one might expect that it is related to packing constraints on the large 
spheres. Other theoretical studies of random sphere parking [23l 121] have calculated the 
maximum number of large spheres iV max that can fit around a small sphere at a given a. 
However, this bound cannot by itself explain why the yield of tetramers can reach 100% 
while that of other clusters, such as trimers or hexamers, cannot. At a given a, it tells us 
only why no clusters larger than iV max (a) can form, but it says nothing about the probability 
of forming smaller clusters with different arrangements. 

Therefore we also examine a different bound, one not previously discussed in the context of 
random sphere parking: the "minimum parking" curve N min (a). N min is the smallest number 
of hard spheres that can be positioned on a smaller sphere such that another sphere cannot 
fit. To understand this bound, consider a simple, one-dimensional analogy to car parking on 
a busy city street, where if a space opens up that is large enough to fit a car, it is filled. The 
minimum parking number occurs when all drivers have been equally inconsiderate, leaving 
spaces between their parked cars that are all slightly too small for another car to fit. This 
lower bound is meaningful only at long times, when all available parking spaces have been 
filled. The long-time limit holds also for our experiments and simulations, which we carry 
out until the average cluster size has saturated. 

Whereas the upper bound iV max (a) is straightforwardly related to solutions of the well- 
known spherical packing problem [HUE?], the calculation of the lower bound N min (a) requires 
a different approach. In our clusters, the distance between the centers of any two big spheres 
must be at least 2i?bi g - Consider then a sphere of radius (-R S maii + Rug) that circumscribes 
the centers of the parked spheres. If this sphere is completely covered with N circles of 
radius 2_R big , it will be impossible to add an (N + l) th large sphere. We are led naturally to 
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the spherical covering problem, a problem with a rich history in mathematics. Like spherical 
packings, spherical coverings are solutions to an extremum problem: they are arrangements 
of TV points on a sphere that minimize the largest distance between any location on the 
sphere surface and the closest point [13]. But unlike spherical packings, spherical coverings 
need not correspond to arrangements of non-overlapping spheres. We therefore solve for the 
minimum parking curve by examining the solutions to the spherical covering problem [27] 
at each N and manually verifying that they correspond to non-overlapping configurations 

Our analytical results for the bounds reveal why a c is a special point: it is the only non- 
trivial point where the calculated maximum and minimum parking curves come together 
(Figure [3]). Analytically we find the location of the critical value to be a c = (1 + y/2) ~ 2.41, 
very close to the values where the experimental distributions are peaked. At a slightly larger 
than this value, the minimum parking configuration corresponds to two spheres placed at 
opposite poles (iV min = 2), and the maximum N is obtained by first parking three large 
spheres next to one another, so that there is room for one more sphere to park (iV max = 4). 
At a slightly smaller than a c , the big spheres can park along orthogonal axes about the small 
sphere to make an octahedron (iV max = 6). The minimum N is obtained by placing four 
spheres as far from each other as possible, so as to make the addition of a fifth impossible 
(-^min = 4). Thus as we increase a through a c , iV max goes from 6 to 4 and N min from 4 to 2, 
and the two curves become infinitesimally close. 

The parking process is therefore geometrically constrained to yield clusters with exactly 
N = 4 particles in the limit a — > a c . A simple geometric argument sheds some light on 
this result. At a c there is always room for four large spheres to park. Parking more spheres 
requires that at least three park precisely along a great circle of the smaller particle, but 
the probability of this happening randomly is zero. Thus irreversible binary aggregation, a 
stochastic process, has a deterministic feature at the critical size ratio: although the space 
between the large spheres can vary, all clusters must be tetramers. Our numerical approach 
confirms that the statistical dispersion in the cluster size distribution vanishes at a c , as 
shown in Figure [3} 

The experimental and simulated distributions differ slightly due to two effects. First, the 
measured sizes tend to be smaller than the simulated ones because a few parking spaces 
remain unfilled even at long times. This effect is more pronounced for larger spheres, which 
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FIG. 3. iVmax (solid gray) and iV m i n (black) as functions of a. Cluster images show sphere 
configurations at discontinuities of these curves. Average cluster sizes from simulations (dashed 
gray) and experiments (blue and red data points) are shown. We characterize the statistical 
dispersion in each distribution by the average absolute deviation from the median, indicated by 
dotted light gray lines for simulations and vertical bars for experiments. 

diffuse more slowly and encounter the small spheres less frequently. The systems most 
affected are the electrostatic ones at a = 3.06 and 4.29. Second, the experimental size 
ratios can vary by 5% due to polydispersity. Both of these factors increase the width 
in the experimental distributions and diminish the achievable yield of tetramers near a c . 
Nevertheless, the experimental data indicate that near a c a tetramer yield of at least 90% 
is possible, and, although the model does not account for the finite range of the interactions 
or the diffusivity of the particles, it is useful for predicting cluster size distributions in two 
very different colloidal systems. 

These results have both fundamental and practical consequences. On the fundamental 
side, the particle size ratio could affect the jamming threshold in bulk packings of bidisperse 
spheres. Previous simulations of these systems have shown that the distribution of coordina- 
tion numbers also depends on the size ratio [28] and may be modeled using random parking 
[25] . This contrasts with dense atomic systems like metallic glasses [T5J HH] in which the 
atoms have some freedom to rearrange locally. In these systems packing constraints may 
explain structure and coordination better than parking arguments. 

On the practical side, this random aggregation process is a simple way to mass produce 
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tetrahedral clusters in theoretically 100% yield. Although the tetrahedra we produce are 
irregular in that the distance between the large spheres can vary, it may well be possible to 
form large quantities of symmetric tetrahedra simply by shrinking the small spheres after 
the tetramers have formed [29]. Furthermore, although the yield will approach 100% only 
for dimers and tetramers, the yield of any iV-particle cluster can be maximized by choosing 
the appropriate size ratio. For instance, the yield of octahedral clusters, also promising 
candidates for building metamaterials [10J, may surpass 70% at a = 1.42. 

The size ratio in binary colloidal systems thus emerges as a valuable control parameter 
for directed self-assembly. Moreover, because it does not require precise control over the 
interactions, random parking offers a robust and simple way to make colloidal clusters that 
are more monodisperse than those prepared through other methods [2J. 
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ELECTROSTATIC INTERACTION EXPERIMENTS 



Charged colloidal polystyrene spheres were purchased from Invitrogen as "IDC surfactant- 
free latex" in batches as listed in Table El 

Mean diameter Surface functionality Fluorescent? Surface charge 

0.49 pm carboxylate-modified latex (CML) yes —262 pC/cm 2 
0.95 |mi amidine no +23.7 pC/cm 2 

1.2 pm amidine no +18.2 pC/cm 2 

1.5 um aldehyde-amidine no +18.2 uC/cm 2 

2.1 um amidine no +30.2 uC/cm 2 

TABLE II. Charged particles used in electrostatic system experiments. Values for surface charge 
from data sheets provided by manufacturers. 

A 100 pL sample of each colloid was diluted to 1% weight by volume. This was then 
vortexed for a few seconds and bath-sonicated for 10 seconds. We cleaned the particles by 
centrifuging and redispersing them several times in deionized (DI) water, using the following 
wash procedure. 

1. Colloids were centrifuged for 5 minutes at 6600g. 

2. Supernatant was removed and 190 pL of DI water were added to each sample. 

3. Samples were vortexed for 5 seconds each and then bath-sonicated for 10 seconds. 

We performed six wash cycles. After the last centrifugation, the supernatant was replaced 
with 40 pL of DI water, rather than 190 pL as before. Then 50 pL of 20 mM NaCl were 
added to each sample to achieve an overall salt concentration of 10 mM. This screens the 
repulsion between like particles before mixing. 

We prepared mixtures of the positively and negatively charged particles such that each 
mixture contained one batch of positively charged particles at 1% w/v. In each mixture, 
the number ratio of the large (positively charged) to small (negatively charged) spheres was 
100 : 1. The salt concentration in each mixture was 10 mM NaCl. Each mixture consisted 



of large and small particles with a different size ratio a, as listed in Table [HI 
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a Large particles 



Small particles 



1.94 0.95 urn amidine (+) 0.49 \im CML (-) 

2.45 1.2 \im amidine (+) 0.49 urn CML (-) 

3.06 1.5 um aldehyde-amidine (+) 0.49 um CML (-) 
4.29 2.1 \im amidine (+) 0.49 um CML (-) 

TABLE III. Size ratios and components of binary mixtures of charged colloids. 

Each mixture was stored in a micro-centrifuge tube and vortexed at 3,000 RPM, bath- 
sonicated for 20 seconds, and then mounted on a Glas-Col Rugged Rotator to tumble slowly 
at 4°C. Each mixture tumbled for at least three days before observation to reduce the effects 
of sedimentation on local particle concentrations throughout the mixture. To make it easier 
to identify and characterize one cluster at a time, we diluted samples of the mixtures to 
0.1% w/v just prior to observing the distribution of cluster sizes. 



Electrostatics control experiment 



In the experiments outlined above, each mixture contained particles with surface charges 
of opposite sign. In a separate control experiment, we mixed particles of two different sizes 
but with surface charges of the same sign. Both components in our control mixture were 



carboxylate-modified latex (CML) colloids with a size ratio a = 2.24, as listed in Table IV 

Mean diameter Surface functionality Fluorescent? Surface charge 

0.49 urn CML yes -262 uC/cm 2 

1.1 urn CML no -31.5 uC/cm 2 

TABLE IV. Colloids used in electrostatic system control experiment. 



These colloids were washed using the procedure outlined above and then mixed in a 
100 : 1 number ratio. After tumbling at 4°C for several days, the cluster size distribution was 
measured. As shown in Figure |4j fewer than 1% of the small particles bind to large particles 
when they have surface charges of the same sign. Nonspecific aggregation is therefore rare 
in the charged colloidal systems. 
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CML control: diameter ratio 2.24 




Number of large spheres in cluster, N 



FIG. 4. Cluster size distribution in a mixture of 1.1 u,m non-fluorescent CML particles and 
0.49 u,m fluorescent CML particles in a 100 : 1 number ratio, showing that particles with surface 
charge of the same sign rarely form clusters. 

Experiments without salt 

In another set of experiments, mixtures of charged colloidal particles were prepared as 
described above, but without salt. The number ratio of positively charged spheres to nega- 
tively charged spheres was again 100 : 1, and several size ratios a were investigated, as listed 
in Table |Vj As in the experiments with 10 mM NaCl, we tumbled the mixtures for several 
days before measuring the distribution of clusters. 

a Large particles Small particles 

1.09 1.2 urn amidine (+) 1.1 u,m CML (-) 

1.36 1.5 ixm aldehyde-amidine (+) 1.1 u,m CML (-) 
1.94 0.95 urn amidine (+) 0.49 urn CML (-) 

2.45 1.2 urn amidine (+) 0.49 urn CML (-) 

3.06 1.5 um aldehyde-amidine (+) 0.49 um CML (-) 
4.29 2.1 urn amidine (+) 0.49 um CML (-) 

TABLE V. Size ratios and components of binary mixtures without salt. 
Figure [5] shows that the average cluster sizes in these mixtures were smaller than the 
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average sizes predicted from simulation and those observed in mixtures containing 10 mM 
NaCl. For instance, at a = 2.45 with 10 mM NaCl the average cluster size is N = 3.9, 
but when there is no salt in the system the average cluster size is N = 2.7. For the four 
size ratios for which there is data both without salt and with 10 mM NaCl, we found that 
clusters are 20% to 35% smaller when no salt is added. 
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FIG. 5. Average cluster sizes from simulations (dashed dark gray) and electrostatic experiments 
with (red data points) and without (cyan) salt. Widths of the cluster size distributions are indicated 
by dotted light gray lines for simulations and vertical error bars for experiments. 



These experiments show that electrostatic repulsion affects the cluster assembly. This 
observation is consistent with other recent experiments [30J showing that the cluster size 
distribution in binary mixtures depends on ionic strength when the salt concentration is less 
than 10 mM. In our mixtures with 10 mM NaCl, the Debye length is approximately 3 nm, 
very small compared to the particle sizes, so the large spheres do not interact with each 
other except at small distances. The random parking model should be more appropriate for 
systems like these where the interaction range is much smaller than the particle size. This is 
because the random parking model assumes no interactions between the particles except for 
a hard-core repulsion and irreversible binding on contact between spheres of two different 
types. 
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DNA-COLLOID EXPERIMENTS 



In another set of experiments, we mixed small and large spheres labeled with comple- 
mentary 65-base ssDNA oligonucleotides purchased from Integrated DNA Technologies: 

• Sequence A: 5'-biotin-51xT-TGTTGTTAGGTTTA-3' 

• Sequence B: 5'-biotin-51xT-TAAACCTAACAACA-3' 

The oligonucleotides terminate with a biotin group, which allows us to graft them to 
streptavidin-coated polystyrene particles using a protocol from Dreyfus et al. [2T] . The 
streptavidin-coated polystyrene particles are purchased from Bangs Laboratories, Inc. with 
the following diameters: 

• 0.21 [xm, fluorescent, to be coated with sequence B 

• 0.39 |xm, fluorescent, to be coated with sequence B 

• 0.51 |xm, fluorescent, to be coated with sequence B 

• 0.97 |xm, non- fluorescent, to be coated with sequence A 

DNA strands were dissolved in water in 20 \xM. concentration. 20 \xL of this DNA solution 
were mixed with 10 \\L of 1% w/w streptavidin-coated polystyrene particles with 120 \xL of 
phosphate buffer in a 1.7 mL propylene micro-centrifuge tube. Each 50 mL batch of buffer 
contained 0.0128 g KH 2 P0 4 , 0.0707 g K 2 HP0 4 , 0.1467 g NaCl, and 0.250 g F108 surfactant 
in 50 mL of deionized water. The buffer was filtered through a 0.2 |j.m membrane before use. 
A separate batch containing no salt was also prepared so that the salt concentration could 
be adjusted by combining the two. Mixtures were vortexed at 3,000 RPM for 5 seconds and 
bath-sonicated for 10 seconds. They were incubated at room temperature for 30 minutes to 
allow the ssDNA to graft to the surface of the particles. Then we washed the colloids using 
the following procedure: 

1. Colloids were centrifuged for 3 minutes at 12,000g. 

2. Supernatant was removed and 100 |oL of 50 mM NaCl buffer were added to each 
sample. 
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3. Samples were vortexed for 5 seconds each at 3,000 RPM and then bath-sonicated for 
10 seconds. 



We washed each colloid three times to remove excess DNA from the system, and then 
incubated each at 55°C for 30 minutes. We then washed three more times, incubated at 55°C 
for another 30 minutes, and washed three times again. At this point the salt concentration 
in the buffer was adjusted to 20 mM. The A- and B-labeled particles were mixed in a 100 : 1 
(large : small) number ratio such that the larger particles were at a volume fraction of about 



0.1%. Three separate mixtures were prepared, as listed in Table VI 



a Large particles Small particles 

1.90 0.97 pm, sequence A 0.51 pm, sequence B 
2.49 0.97 pm, sequence A 0.39 pm, sequence B 
4.62 0.97 pm, sequence A 0.21 pm, sequence B 

TABLE VI. Size ratios and components of mixtures with DNA-driven interactions. 

After the mixtures had been prepared, we followed the same procedure that we used for 
the charged colloid system. 



DNA-colloid control experiment 

In the experiments outlined above, each mixture contained particles labeled with comple- 
mentary DNA strands. In a separate control experiment, we mixed particles of two different 
sizes but labeled with identical ssDNA. 

Both components in our control mixture were streptavidin-coated polystyrene colloids 
labeled with sequence A. We used 0.97 pm (non- fluorescent) and 0.51 pm (fluorescent) par- 
ticles, yielding a size ratio a = 1.90. The particles were functionalized with DNA using the 
same procedure described above and then mixed in a 100 : 1 number ratio. After tumbling at 
4°C for several days, the cluster size distribution was measured. Figure [6] shows that fewer 
than 2% of the small spheres bind to large spheres when they are coated with the same 
DNA sequence. The low amount of nonspecific aggregation is expected, since we designed 
sequence A to have a negligible amount of self-hybridization even at 0°C 
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FIG. 6. Cluster size distribution in a mixture of 0.97 u,m non-fluorescent particles and 0.51 u,m 
fluorescent particles in a 100 : 1 number ratio, coated with the same non-self-complementary DNA 
sequence (A). 



MEASUREMENT OF CLUSTER SIZE DISTRIBUTION 

To measure the distribution of cluster sizes of a particular mixture, we placed a 5 \xL, 
0.1% w/v sample between two cover slips and then sealed it at the edges with UV-curable 
epoxy (Norland Optical Adhesive 61). We observed each sample under differential interfer- 
ence contrast with a 100X oil-immersion objective on a Nikon Eclipse TE2000-E inverted 
microscope. We used fluorescence to identify the small spheres. 

The distribution was obtained by counting the number of clusters at each N. N = is a 
small sphere without any large spheres adsorbed, N = 1 a. small sphere attached to a single 
large sphere, etc. We also counted clusters containing multiple small spheres and clusters 
formed through non-specific aggregation. The number of such clusters is small compared to 
the total number of clusters counted. Thus, we do not include these in the histograms. 

Our counting procedure is designed to avoid errors from double counting. We begin by 
counting clusters in the field of view (FOV) on the microscope, scanning through z to find 
clusters that may initially be out of focus. The Brownian motion of each cluster eventually 
brings all particles into view. We can therefore determine the number of large and small 
spheres in each cluster through direct observation. Once we have recorded all clusters in 
the FOV, we translate the FOV to another part of the sample located more than one FOV 
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away. We repeat this process to build up the histogram, rastering the FOV over the sample. 
We count more than 60 clusters at each size ratio and more than 120 clusters at size ratios 
greater than 2. The entire histogram for a given mixture is recorded in one session on the 
microscope. 



ALGORITHM FOR SIMULATING RANDOM SPHERE PARKING 

We simulate the cluster distribution using an algorithm based on Monte Carlo trials. 
These consist of two stages: 

Coarse stage We repeatedly try to insert a disc of radius r = i?bi g / (-Rbig + -Rsmaii) on the 
surface of the unit sphere. The center of the disc is randomly and uniformly distributed on 
the sphere. If the disc overlaps any discs that are already "parked", we reject it; otherwise 
we add it to the list of parked discs. After a fixed number of consecutive rejections, (typically 
^coarse = 10 4 ), we switch to the fine stage. 

Fine stage We compute the remaining regions of possible insertion, and then try to 
insert an additional disc. We choose the center of the disc randomly and uniformly from 
these regions. 

To find the regions of possible insertion, we first compute the arcs that form their bound- 
aries: Given a central disc of radius r, neighboring disc centers must lie on or outside a 
concentric circle of radius 2r. We erase arcs that lie in the interior of the concentric circles 
of radius 2r about each neighboring disc. If the central circle contains any unerased arcs, 
these are added to our roster of arcs. 

We do this for each parked disc center, and then stitch together the remaining arcs to 
form the regions to which another disc center could be added. To find the area of each 
region, we inscribe the region in a circumcircle and use a Monte-Carlo integration method 
to determine the ratio of the area of the region to the area of its circumcircle. 

The algorithm terminates when the list of remaining arcs in the fine stage is empty, at 
which point the final number of inserted discs is recorded. 
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CALCULATING BOUNDS ON CLUSTER SIZE DISTRIBUTION 

The upper and lower bounds N max and iV min on the cluster size distribution can be calcu- 
lated from spherical codes [27J corresponding to known solutions to the problems of spherical 
packings (N max ) and spherical coverings (iV m i n ). Spherical packings are arrangements of N 
points on a unit sphere that maximize the smallest distance between any two of them [131 E7J • 
For a given a we determine N max (a) by looking up the spherical packing [27J with the largest 
N for which the minimal distance between points is at least 2a/(l + a). Our calculation of 
-^mm(«) is more involved, as explained below. 



Connection between spherical covering and minimum parking 

Here we verify that the best known coverings are also optimal solutions to the minimum 
parking problem. For a given configuration of n "parked" points on a unit sphere, let the 
covering radius be the maximum distance between any point on the sphere to the nearest 
parked point. Let the packing radius be (half) the minimum of the pairwise distances 
between the parked points. If the parked points represent centers of circles with some radius 
r, it is impossible to add another circle without overlapping a parked one if and only if the 
covering radius is less than 2r. The circles that are already parked do not overlap provided 
that the packing radius is greater than r. Therefore, if we are given an optimal solution to 
the covering problem, that is, one that minimizes the covering radius for a given n, it will 
also be the optimal minimal parking configuration if the covering radius is less than 2 times 
the packing radius. 

We manually verify that the best known coverings satisfy this constraint for n — 4, . . . , 130 
by calculating the packing radius for each optimal known covering, obtained from ref. 



Table VII shows the covering radius, packing radius, and difference (in degrees). The dif- 
ference is always positive, confirming our statement. 

n 2 * packing radius covering radius 2* (packing ra- 
dius) - (covering 
radius) 



4 109.47 70.53 38.94 

5 90.00 63.43 26.57 
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6 


90.00 


54.74 


35.26 


7 


72.00 


51.03 


20.97 


8 


61.76 


48.14 


13.62 


9 


68.97 


45.88 


23.09 


10 


65.53 


42.31 


23.22 


11 


50.65 


41.43 


9.22 


12 


63.43 


37.38 


26.06 


13 


46.23 


37.07 


9.16 


14 


52.58 


34.94 


17.64 


15 


45.67 


34.04 


11.63 


16 


50.48 


32.90 


17.58 


17 


41.63 


32.09 


9.54 


18 


45.53 


31.01 


14.51 


19 


40.73 


30.37 


10.36 


20 


40.01 


29.62 


10.39 


21 


39.45 


28.82 


10.62 


22 


40.70 


27.81 


12.89 


23 


38.99 


27.48 


11.51 


24 


36.67 


26.81 


9.86 


25 


36.75 


26.33 


10.42 


26 


35.12 


25.84 


9.27 


27 


38.06 


25.25 


12.81 


28 


35.87 


24.66 


11.21 


29 


33.86 


24.37 


9.50 


30 


31.18 


23.88 


7.30 


31 


29.94 


23.61 


6.33 


32 


37.38 


22.69 


14.69 


33 


26.61 


22.59 


4.02 


34 


30.82 


22.33 


8.49 


35 


27.98 


22.07 


5.90 



36 


28.78 


37 


31.22 


38 


30.31 


39 


30.73 


40 


30.13 


41 


27.71 


42 


28.34 


43 


27.27 


44 


26.36 


45 


25.15 


46 


29.11 


47 


25.38 


48 


27.70 


49 


24.93 


50 


28.01 


51 


24.30 


52 


24.24 


53 


22.56 


54 


25.58 


55 


24.06 


56 


24.89 


57 


23.98 


58 


23.72 


59 


23.69 


60 


24.67 


61 


20.50 


62 


20.43 


63 


21.60 


64 


21.34 


65 


20.78 



21.70 


7.08 


21.31 


9.91 


21.07 


9.24 


20.85 


9.87 


20.47 


9.66 


20.32 


7.39 


20.05 


8.29 


19.84 


7.43 


19.64 


6.72 


19.42 


5.73 


19.16 


9.96 


18.99 


6.38 


18.69 


9.01 


18.59 


6.33 


18.30 


9.71 


18.20 


6.10 


18.05 


6.19 


17.88 


4.68 


17.68 


7.90 


17.52 


6.54 


17.35 


7.54 


17.18 


6.80 


17.02 


6.70 


16.90 


6.79 


16.77 


7.90 


16.64 


3.86 


16.49 


3.94 


16.37 


5.23 


16.19 


5.15 


16.11 


4.66 



66 


21.61 


67 


21.32 


68 


21.74 


69 


20.41 


70 


21.47 


71 


21.46 


72 


23.06 


73 


17.51 


74 


17.95 


75 


20.02 


76 


19.08 


77 


22.03 


78 


19.14 


79 


19.77 


80 


19.94 


81 


19.07 


82 


20.40 


83 


17.08 


84 


19.93 


85 


19.92 


86 


19.89 


87 


18.70 


88 


19.80 


89 


19.76 


90 


19.83 


91 


17.18 


92 


17.83 


93 


18.75 


94 


19.15 


95 


18.54 



15.96 


5.66 


15.86 


5.46 


15.72 


6.02 


15.60 


4.82 


15.50 


5.98 


15.39 


6.07 


15.14 


7.92 


15.12 


2.39 


15.03 


2.92 


14.95 


5.07 


14.85 


4.23 


14.74 


7.29 


14.66 


4.49 


14.56 


5.21 


14.45 


5.49 


14.38 


4.69 


14.29 


6.11 


14.22 


2.86 


14.12 


5.81 


14.05 


5.88 


13.96 


5.93 


13.88 


4.81 


13.79 


6.01 


13.71 


6.05 


13.62 


6.21 


13.56 


3.62 


13.49 


4.34 


13.43 


5.32 


13.35 


5.81 


13.29 


5.25 



96 


19.21 


13.21 


6.00 


97 


18.38 


13.14 


5.24 


98 


18.99 


13.06 


5.93 


99 


19.06 


13.00 


6.06 


100 


18.58 


12.94 


5.64 


101 


18.67 


12.87 


5.80 


102 


17.63 


12.81 


4.82 


103 


17.51 


12.74 


4.77 


104 


18.49 


12.67 


5.82 


105 


18.94 


12.62 


6.32 


106 


17.66 


12.56 


5.10 


107 


18.62 


12.50 


6.13 


108 


18.10 


12.43 


5.67 


109 


18.23 


12.38 


5.85 


110 


18.40 


12.30 


6.10 


111 


18.70 


12.25 


6.45 


112 


18.71 


12.19 


6.52 


113 


18.28 


12.15 


6.14 


114 


16.87 


12.10 


4.78 


115 


17.42 


12.05 


5.37 


116 


16.27 


11.99 


4.28 


117 


17.67 


11.94 


5.73 


118 


17.71 


11.89 


5.83 


119 


17.09 


11.84 


5.25 


120 


17.25 


11.79 


5.47 


121 


17.56 


11.73 


5.83 


122 


17.70 


11.68 


6.02 


123 


14.38 


11.64 


2.74 


124 


14.14 


11.59 


2.55 


125 


17.38 


11.54 


5.84 



126 17.07 11.49 5.58 

127 16.54 11.45 5.09 

128 16.50 11.41 5.09 

129 16.98 11.36 5.62 

130 16.95 11.32 5.63 

TABLE VII: Packing and covering radii for different clus- 
ter sizes n, measured in degrees, and the differences be- 
tween them. Because the differences are always positive, 
we conclude that the optimal covering configuration is also 
the optimal minimal parking configuration, at least for N G 
{4,..., 130}. 
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